Vortex lattice structures and pairing symmetry in Sr2Ru04 
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Recent experimental results indicate that superconductivity in Sr2Ru04 is described by the p-wave repre- 
sentation of the D4h point group. Results on the vortex lattice structures for this representation are presented. 
The theoretical results are compared with experiment. 



1. Introduction 

The oxide Sr2Ru04 has a structure similar to 
high Tc materials and was observed to be su- 
perconducting by Maeno et al. in 1994 ||l[. R 
has been established that this superconductor is 
not a conventional s-wave superconductor: NQR 
measurements show no indication of a Hebel- 
Slichter peak in 1/TiT |], and is strongly 
suppressed below the maximum value of 1.5 K 
by non-magnetic impurities More recent ex- 
periments indicate an odd parity gap function of 
the form d(k) — z{'qxkx + Vyky). The Knight 
shift measurements of Ishida et al. Q reveal that 
the spin susceptibility is unchanged upon enter- 
ing the superconducting state; this is consistent 
with p-wave superconductivity (as predicted by 
Rice and Sigrist ||^). Furthermore, these mea- 
surements were conducted with the applied field 
in the basal plane. Since the orientation of the 
gap function (that is d) is orthogonal to the 
spin projection of the Cooper pair j6j, these mea- 
surements are consistent with the gap function 
aligned along the z direction. The /xSR experi- 



ments of Luke et. al. have revealed spontaneous 
fields in the Meissner state 0. This indicates 
that the superconducting order parameter must 
have more than one component 1^,0]. Naively, 
the observation of spontaneous fields leads to the 
conclusion that the superconducting gap function 
breaks time reversal symmetry (T) and therefore 
must have the form d(k) — z(kx ±iky). However, 
the muons probe inhomogeneities in the super- 
conducting state since any bulk magnetic fields 
are screened in the Meissner state. Consequently, 
T need only be broken in the vicinity of the inho- 
mogeneities. For example, it has been pointed out 
in Ref. that a gap function given by rf = (1, ±1) 
can give rise to domain walls that break time re- 
versal symmetry. In principle this can also lead 
to a /xSR signal as seen by Luke et al.. Conse- 
quently, these two experiments lead to the gap 
function d(k) = z{'qxkx -'r'nyky), with the specific 
form of the order parameter rix and rjy undeter- 
mined. In this talk, the vortex lattice structures 
arising in the phcnomenological theory of this Eu 
representation are examined. First, the theory 
is examined for an applied magnetic field in the 
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basal plane along one of the two-fold symmetry 
axes. It is shown that multiple vortex lattice 
phases are generic to this representation. Then 
the weak-coupling limit of the theory is examined 
for the applied magnetic field along the c-axis. It 
is shown that a square vortex lattice results for 
parameters relevant to Sr2Ru04; consistent with 
the experimental observation of a square vortex 
lattice . The theory of the square vortex lattice 
is compared with the observed results and good 
agreement is found. 

2. Free Energy 

The dimensionless free energy for the £"„ rep- 
resentation of D^h is given by ^ 

/= -\ri\^ + \ri\y2 + P2iv.V*y-VyV*.)y'2 
+fi3\r].\^\%\^ + \D.V.\^ + \Dy7Jy\^ 

+ ^i2{\DyV.\^ + \D,lly\^) (1) 

+H3[{Dxrix){Dyriy)* + h.c] 

+ K4iDyVx)iDxVy)* + h.c] + / (Stt) . 

where Dj = Vj - ff^j, h = V x A, and 
A is the vector potential. The stable homo- 
geneous solutions are easily determined. There 
are three phases: (a) — (l,i)/-\/2 (/32 > 
and /32 > /33/2), (b) 77 = (1,0) (/Jg > and 
P2 < (33/2), and (c) rf = (l,l)/x/2 (P^ < and 
P2 < 0). Phase (a) is nodeless and phases (b) 
and (c) have line nodes. For some of the dis- 
cussion in this paper the Ginzburg Landau co- 
efficients are determined within a weak-coupling 
approximation in the clean limit. In this case, 
taking for the Eu REP the gap function de- 
scribed by the pseudo-spin-pairing gap matrix: 
A = ilviVx/ y/ {vl)+r]2i>y/\/ (w^)]^^^^, where the 
brackets () denote an average over the Fermi sur- 
face and ai are the Pauli matrices, it is found 
that [|| P2 — H2 — H3 = Hi — {i' + — v) and 
/?3 = 4i^/(3 — v) where 



Note that < 1 and v — ioY a. spherical or 
a cylindrical Fermi surface. Also note that fj oc 
(1, i) is the stable homogeneous state for all v. 



3. Magnetic field in the basal plane 

For the Eu model, symmetry arguments im- 
ply that the vortex lattice phase diagram con- 
tains at least two vortex lattice phases for mag- 
netic fields applied along at least two of the four 
two-fold symmetry axes |9j. To demonstrate this, 
consider the magnetic field along the x direction 
{x is chosen to be along the crystal a axis) and 
a homogeneous zero-field state 77 cx (1,?). The 
presence of a magnetic field along the x direction 
breaks the degeneracy of the (77^;, rjy) components, 
so that only one of these two components will or- 
der at the upper critical field [e.g. fJ oc (0,1)]. 
It has been shown for type II superconductors 
with a single component order parameter that 
the solution is independent of x [|lo| so that ax 
(a reflection about the x direction) is a symmetry 
operation of the ff cx (0, 1) vortex phase. Now 
consider the zero-field phase ff oc (1,?); cr^ trans- 
forms (1, i) to (—1, i) e*'^(l, i) where tp is phase 
factor. This implies that is not a symmetry 
operator of the zero-field phase. It follows that 
there must exist a second transition in the fi- 
nite field phase at which rjx becomes non-zero. 
Similar arguments hold for the field along any of 
the other three two-fold symmetry directions in 
the basal plane. Consequently a zero-field state 
ff cx (l,i) must exhibit two vortex lattice phases 
when the field is applied along any of the four 
two-fold symmetry axes. Similar arguments for 
ff cx (1,0) or 77 cx (1,1). imply that there must 
exist at least two vortex lattice phases for only 
two of the four two-fold symmetry axes. For a 
zero-field state 77 = (1,0) fields along the (1,1) or 
the (1,-1) directions will result in two vortex lat- 
tice phases (for the (1, 0) and the (0, 1) directions 
multiple vortex lattice phases may exist, but are 
not required by symmetry). For a zero-field state 
rf cx (1, 1), two vortex lattice phases will exist for 
fields along the (1,0) or the (0, 1) directions. 

An analysis of the free energy of Eq. Q re- 
veals additional information about the phase di- 
agram HJllI]. Assuming the large k limit (note 
K « 30 for this field orientation in Sr2Ru04), 
the vector potential can be taken to be A = 
_ff2;(sin 0, — cos (/>, 0) {(j) is the angle in the basal 
plane that the applied magnetic field makes with 
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the X direction). The component of D along the 
field is set to zero. The upper critical field found 
by this method exhibits a four-fold anisotropy in 
the basal plane ||;^. For concreteness, the field 
is taken along the x direction {(j) = 0). Intro- 
ducing the raising and lowering operators n± = 
q{\J K2I K.'^Dy ± iD^)/y/2 with = ^/k^Jk^/H, 
the gradient portion of the free energy can be 
written as 

fgrad= ^/T^H {7^*[l + 2N]r^^ + T^;[{\ + 

{l + 2N) + {^^-\){T\l+TP)]r^y}. 

Assuming that K2 < 1 the first transition is given 
by the standard hexagonal Abrikosov vortex lat- 
tice solution for r^x- To find the second transition 
the lowest eigenstate oi{\ + ^){l + 2N) + {^- 

1) (e2'^n2_-he-2*^n?_) must be found (where 9 has 
been introduced so that the vortex lattice can be 
rotated with respect to the ionic lattice) and a 
vortex lattice solution for rjy must be constructed 
from this. Note that the zeroes of the rjx lattice 
and the zeroes of the 7]y lattice need not coin- 
cide ijl^. This method yields for the ratio of the 
second transition {H2) to the upper critical field 

{Hc2) 

H2 _ Pa — Pm 

where (3a = 1.1 596, = (1 - /^a + fiz)Si - 

l/?2g2|, ^1 ^ \m?\m?l{\m? \m?), = 

{'ni'H'iY I ilViV hap), the over-bar denotes a spa- 
tial average, and /3m must be minimized with re- 
spect to 9 and the displacement between the zeros 
of the rjx and the rjy lattices. Diagrams of the pre- 
dicted vortex structures are shown in Fig. ^ and 
the phase diagram for the weak-coupling limit is 
shown in Fig. ^ Two vortex lattice configura- 
tions are found to be stable as a function of v. 
For —0.23 < V < 1 (Phase 1), the displacement 
between the the zeroes of the r]x and r]y lattices is 
one half of a hexagonal lattice basis vector and 9 
varies with v (and in general with field and tem- 
perature) so that the vortex lattice is not aligned 
with the ionic lattice. For — l<i^< —0.23 (phase 

2) , the r]x and rjy lattices coincide and a vortex 
lattice basis vector lies along the i/ direction. Kita 
|Tl| has argued that for phase 1, at some field be- 
low H2 there will be a first order transition from 
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Figure 1. Stable vortex lattice structures found 
for the field along x. The open (closed) circles 
correspond to zeroes of the rjx {%) lattice. 



phase 1 to phase 2. For the field along x zLy the 
phase diagram is given by replacing ly with —ly. 
Recent experiments of Mao et al. |l^ show some 
support for the theory presented here. 

4. Magnetic field along the c-axis 

The free energy in Eq. |^ in the weak coupling 
limit has been used to determine the vortex lat- 
tice structure for the field along the c-axis for 
fields near the upper critical field |lj]. The main 
conclusion of this analysis is that vortex lattice 
is square for > 0.0114 and that the vortex 
lattice is oriented along the crystal lattice for 

< —0.0114 and rotated 45 degrees from the 
crystal lattice for 1/ > 0.0114. The analysis near 
Hci has also been carried out with the result that 
the vortex lattice at Hd will be hexagonal and 
with increasing field the lattice will continuously 
distort until a square vortex lattice is formed 
p^ ]. A square vortex lattice has been observed 
by Small Angle Neutron scattering (SANS) and 
the orientation implies v > 0.0114 g]. One no- 
table feature of the SANS measurements is the 
observation of Bragg peaks beyond the first Bragg 
peak. The analysis of Eq. |l| has been extended to 
fields below Hc2 and the lowest Bragg peaks of the 
field distribution have been calculated and com- 
pared to experiment. One comparison is shown 
in Fig. H (note that in the SANS measurements 
there is flux pinning which is not included in the 
calculations). The large size of the higher order 
Bragg peaks relative to that expected for a single 



Figure 2. The ratio of the two transition fields 
and 9 (the angle relating the orientation of the 
vortex lattice to the ionic lattice) as a function of 
v. 
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complex order parameter component (Abrikosov) 
theory and the agreement between the theory 
and the experiment gives support to the the- 
ory. However, note that it is possible that an 
Abrikosov theory with sufficiently large non-local 
corrections may also account for the observed re- 
sults. 
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